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Abstract 


The  interaction  between  vorticity  and  a  free  sur¬ 
face,  modeled  as  a  shear  free  boundary  is  studied  us¬ 
ing  a  direct  numerical  simulation  of  an  open  channel 
flow  at  low  Reynolds  number  (Jle7T-®-f//i^/i''=-2340 
wberex/i-is-the  channel'depth-and  (//T-is  the  mean  free-Cu. 
i^urface-velocity)!  As  a  result  of  the  shear  free  nature 
of  the  top  boundary,  only  normal  vorticity  may  ter¬ 
minate  on  it.  The  vorticity  components  parallel  to 
the  top  boundary  must  go  to  zero  as  the  boundary  is 
approached.  The  time  averaged  fluctuating  enstrophy 
balance  equations,  which  are  an  indicator  of  the  level 
of  activity'  of  the  vorticity  field  are  evaluated.  Near 
the  free  surface  the  rate  of  production  and  destruction 
of  enstrophy  is  set  by  the  stretching  and  rotation  of 
fluctuating  vorticity  by  the  fluctuating  velocity  field. 
The  transport  of  enstrophy  by  the  fluctuating  veloc¬ 
ity  has  a  minor  role.  Of  the  the  production  terms, 
the  three  involving  the  stretching  of  vorticity  along 
its  axis  are  the  most  important,  and  a  simple  model 
explaining  this  is  presented. 


■J-  JJ 


u>-  , 


Nomenclature 


U,{X2) 

Uf.  =  m{x2  =  0.0) 

U*  =  s/TuialllP 
/•  =  u/u* 

<•  = 


Instantaneous  velocity 
Fluctuating  velocity 
Instantaneous  vorticity 
Fluctuating  vorticity 
RMS  fluctuating  velocity 

RMS  fluctuating  vorticity 

Mean  velocity 
Mean  free  surface  velocity 
Friction  velocity 
Viscous  lengthscale 
Viscous  timescale 
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uf  =•  u,/u‘ 

C  =  un- 

t*  =  t/f 
(u(x))  — 

h  =  Lr, 

Lxj 


P 

Fn^  Vj./s/g^ 
Re  =  Uh/u 
Re-  =  u‘h/u 
Ref,  =  Uj,hlu 


Non-dimensional  velocity 
Non-dimensional  l**ngth 
Non-dimensional  tim>’ 
o(x)  Averaging  operator 
Channel  height 
Stre?r.iwise  box  length 
Transverse  box  length 
Kinematic  viscosity 
density 

Froude  number 
Reynolds  number 
Wall  Reynolds  number 
Free  surface  Reynolds  no. 


1.  Introduction 

The  total  instantaneous  enstrophy  of  a  fluid  is 
defined  as  the  square  of  the  instantaneous  vorticity. 

Enstrophy  can  be  determined  as  well  from 
the  mean  vorticity.  and  the  fluctuating  vortu-- 

ity,  This  paper  will  be  limited  to  an  anal>>ii' 

of  the  fluctuating  vorticity.  The  dominate  motiva¬ 
tion  for  studying  the  enstrophy  balance  equations  is 
the  natural  interpretation  of  enstrophy  as  a  measur** 
of  vortical  activity.  The  no-slip  boundary,  where  thi* 
majority  of  the  fluctuating  vorticity  is  produced  is  th** 
source  of  the  turbulent  kinetic  energy.  .Much  of  i!;- 
turbulent  kinetic  energy  is  believed  to  be  associat*- 1 
with  coherent  structures  and  a  quasi-regular  turbu 
lent  burst  cycle.  If  it  is  assumed  that  coh>-r-  ui  voru-' 
ity  is  an  important  feature  of  these  structures,  th-n 
enstrophy  is  a  good  tool  for  analysis.  Although  a  ch-  ir 
connection  between  vorticity  and  coherent  siructuf' 
or  events  has  yet  to  be  developed,  the  concept  of  th** 
vortex  loop  has  often  invoked  to  explain  the  visu.d 
ization  of  the  results  of  experiments  or  direct  numer¬ 
ical  simulations  (DNS).  The  connection  between  vor¬ 
ticity  and  turbulent  structures  may  be  equally  vah  1 
at  the  free-surface.  As  will  be  discussed  below,  the 
conditions  at  the  free-surface  require  that  vorlicitv 
tangential  to  the  surface  vanish  as  the  surface  is  ap 
preached.  Simultaneously,  the  surface-normal  gra-h 
ent  of  the  component  of  vorticity  perpendicular  to  t  h* 
surface  must  vanish.  These  are  reflectional  boun  ; 
ary  conditions,  and  are  significantly  diPferent  than  i!.- 


/)'(  I 


condition  of  rotational  lymmetry  about  the  xi-axis 
(aligned  with  the  channel  centerline)  that  occure  in 
the  closed  channel  case. ' 

Another  important  interpretation  of  enstrophy 
can  be  made  through  its  relation  to  the  isotropic  dis¬ 
sipation  function.  The  isotropic  dissipation  function, 
commonly  used  mk  -  €  turbulence  models,  is  defined 

as; 


dxj  dxj ' 


(1) 


for  homogeneous  turbulence.  However,  the  following 
can  be  shown 


dui  dui  _  I  dui  dui 

dxj  dxj  ~  2  dXj  dxj  2  dz,  dzt 

-  duj  duj  du.  (2) 

~2'dxj  dxi^dxj  dz, 

=  WkW*. 


In  Eq.  2  the  identity  for  homogeneous  incompressible 
turbulence, 


dui  duj  _  d  !  duj  ^  ()  du, 

dxj  dxi  ~  dxj  dxj  dz,  dz,  (3) 

•  =  0, 


has  been  used‘.  Therefore  the  isotropic  dissipation 
function  is  related  to  the  fluctuating  enstrophy  a.s 


(  = 


(4) 


Consequently,  the  study  of  the  enstrophy  balance  can 
be  viewed  as  the  study  of  the  creation  and  destruc¬ 
tion  of  the  isotropic  dissipation  function.  The  enstro¬ 
phy  balance  equation  may  be  a  more  convenient  tool 
to  study  the  dissipation  because  its  terms  have  bet¬ 
ter  physical  interpretations.  Furthermore  the  com¬ 
ponent  equations  UiUi  (no  implied  summation)  can 
be  examined  individually.  As  mentioned  above,  while 
the  isotropic  dissipation  is  not  the  local  dissipation  of 
turbulent  kinetic  energy,  it  is  relevent  for  turbulence 
modeling. 


2.  The  Numerical  Simulation 

The  incompressible  three-dimensional  Navier- 
Stokes  equations  are  solved  for  initial  and  boundary 
conditions  approximating  a  turbulent  open-channel 
flow  of  water  at  a  Reynolds  number  of  2340,  where 
Rth  s  Ui,hlv  and  h  is  the  channel  depth  and  f//i 


Wj •  —I  U/j  .  **1/3  . 


is  the  mean  free-surface  velocity.  Based  on  the  fric¬ 
tion  velocity,  u.  =  the  Reynolds  number  is 

Rt,  =  134.  The  governing  equations  are  recast  into 
a  4'^  order  equation  for  the  vertical  velocity  and  .1 
2"'*  order  equation  for  the  vertical  vorticity  and  the 
continuity  equation  is  solved  explicitly  in  the  recrA 
ery  of  the  streamwise  velocity.  This  formulation  is  l- 
first  discussed  by  Orszag  and  Patera*  and  later 
plementcd  in  a  simpler  form  in  the  large  scale  'lir-  ■ ' 
simulations  by  Kim,  .Moin  and  Moser*.  It  is  ili*.*  i  v' 
ter  formulation  that  will  be  used.  This  method  ui 
volves  the  use  of  the  homogeneous  solutions  of  :!.• 
time  discretized  4‘*  order  equation  to  satisfy  all  'h- 
required  boundary  conditions.  The  equations  are  nu 
merically  solved  after  they  are  Fourier  iransforin-  : 
in  the  streamwise  (xi)  and  spanwise  fz.d  lirecti.iii' 
and  Chebychev  transformed  in  the  vertu-.-xi  .lirecti'  u 
(rj).  The  xj  axis,  scaled  by  the  channel  height  h a.- 
its  origin  at  the  free  surface.  The  channel  bottom 
Z2  =  -1.  The  calculations  are  performed  on  a  grid  ■ ! 
48x65x64  nodes  in  Xi.zj.ra,  respectively,  for  a  n-v- 
lulion  of  26. 3/.  in  the  streamwise  direction  and  13  '21, 
in  the  spanwise  direction  (  /.  =  *^^/^^^)  iThe  t>M.M 

box  size  is  1684/.  x  134/,  x  632/,.  With  the  geometr> 
scaled  by  the  channel  height,  the  height,  length  .an  ! 
width  of  the  channel  are  I,  4>r  and  3t/2,  respect  1  a >  !> 
In  the  free-surface  problem  the  nonlinear  an  -, 
lime-dependent  boundary  conditions  must  be  ■‘aii' 
fled  on  an  unknown  surface  elevation,  t).  The  bou;: : 
ary  conditions  can  be  simplified  considerably  if 
surface  is  not  allowed  to  deflect  (tj  s  0).  The  r  c 
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lid  approximatioa  is  equivalent  to 

tij  =s  0  at  *j  ss  0.  (5) 


To  further  simplify  the  boundary  conditions,  the 
usual  balance  of  normal  stresses^  is  replaced  by 


=  0  at  Z}  =:  0. 


(6) 


Eq.  6  can  be  derived  from  the  zero  tangential  stress 
conditions,  the  continuity  equation,  and  the  rigid  lid 
assumption.  These  conditions  and  the  definition  of 
vorticity  can  also  be  combined  to  form  the  remaining 
free-surface  boundary  condition, 

dui’l 

^=0atz,  =  0,  (7) 

where  the  vorticity  is  non-dimensionalized  by 
t/o  is  the  initial  (<  =  0)  free-surface  velocity. 

At  the  bottom  wall  of  the  channel,  the  no-slip  ve¬ 
locity  boundary  conditions  are  composed.  The  no-slip 
velocity  condition  sustains  the  shear  layer  which  is  re¬ 
sponsible  for  maintaining  the  turbulence.  The  bound¬ 
ary  conditions.are  periodic  on  all  dependent  variables 
in  the  streamwise  and  spanwise  directions. 

After  the  wall  shear  stress  achieved  a  statistically 
steady  state  at  the  correct  value,  and  the  total  stress 
across  the  channel  was  linear,  turbulent  flow  data  was 
acquired.  Forty-two  realizations  of  the  turbulent  flow 
were  saved  during  an  interval  of  approximately  4000t* 
(t*  =  i//u*^).  We  have  compared  single  point  statis¬ 
tics  and  spectra  with  from  the  closed  channel  calcu¬ 
lation  performed  at  the.  Naval  Research  Laboratory'^ 
and,  where  possible,  results  from  the  high  resolution 
NASA  Ames  closed  channel  computation^.  In  gen¬ 
eral,  the  single  point  statistics  in  the  wall  and  buffer 
regions  are  very  similar  to  the  closed  channel  results. 
The  effects  of  the  upper  boundary  are  apparent  in  the 
range  of  0.0  <  2}  <  -0.3  as  compared  to  similar  data 
from  the  closed  channel  calculations. 

The  near  wall  behavior  of  the  rms  vorticity,  wj 
in  Fig.  I  is  similar  to  the  prior  results  of  Kim,  Moin 
and  Moser^.  Near  the  free  surface  the  horizontal  com¬ 
ponents  of  vorticity  vanish,  while  the  vertical  compo¬ 
nent  approaches  a  constant  value  of  =  0.035. 

This  can  be  compared  to  0.067  at  the  centerline  for  a 
low  resolution  closed  channel  simulation  performed  at 
NRL  and  0.042  for  the  results  Kim,  Moin  and  Moser^. 
A  Taylor  series  expansion  about  the  top  boundary  in¬ 
dicates  that  u\  and  both  approach  zero  linearly. 
A  physical  explanation  of  why  I  ^  l>l  ^  I  “  the 
free  surface  is  being  approached  can  be  developed  by 
examining  the  S7^  balance  equations. 


4.  Balance  Eau^^ni 

The  relationship  of  enstropby  to  fluctuating  vor- 
ticity  is  similar  to  the  relationship  of  turbulent  kinetic 
energy  to  fluctuating  velocity.  The  balance  equations 
are  derived  in  the  same  manner  as  the  balance  equa¬ 
tions  for  the  components  of  the  Reynolds  stress  ten¬ 
sor,  The  equations  for  the  mean  and  instantaneous 
vorticity  are  multiplied  the  mean  and  instantaneous 
vorticity,  respectively.  These  equations  are  time  av¬ 
eraged  resulting  in  equations  for  the  mean  enstrophy 
and  the  mean  equation  of  the  instantaneous  enstro¬ 
phy.  Subtracting  the  latter  from  the  former  yields 
the  equation  of  the  mean  fluctuating  enstrophy  .-\s 
in  the  case  of  the  Reynolds  stress  tensor,  the  •“ym- 
tions  for  the  individual  components  of  the  enstrophy 
rTw7,  ^juij ,  Caw?  can  be  formulated.  The  det-ails  of 
the  derivation  of  the  balance  equation  for  the  fluctu¬ 
ating  enstrophy  are  contained  in  Balitit,  Vukosl.iv'-*- 
vie  and  Wallace^.  The  steady  state  balance  equations 
are; 


with  I  s  1,2  or,  3  for  the  component  equations  and 
an  implied  summation  for  the  total  fluctuating  en¬ 
strophy. 

According  toTennekes  and  Lumley*  and  Bahru' 
et  al.,  the  eight  terms  of  the  enstrophy  balance  equ.t- 
tion  have  the  following  interpretation: 

(1)  The  convection  of  the 
fluctuating  enstrophy, 

(2)  Transport  of  fluctuating 
enstrophy  by  fluctuating  velocity, 

(3)  Gradient  production 
of  enstrophy, 

(4)  ^7371^:  production 

by  mean  velocity  gradient, 

(5)  production 

by  fluctuating  velocity  gradient. 

(6)  ;  .Mixed  production, 

(7)  ;  Viscous  diffusion, 

(«)  “Wf:  Viscous  dissipation. 

The  interpretation  of  terrns  1-6  can  b**  -  x 
panded  by  examining  the  non-linear  te^m  of  the  x  - 
ticity  equation  from  which  they  were  derived; 

V  X  u  X  w  s  (  w  •  V)n  -  (u  •  V)  w 
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The  first  term  U  the  smount  of  init&ntaoeoui  stretch¬ 
ing  or  rotation  of  the  vortlcity  by  the  gradient  of  the 
velocity  field.  Following  Batchelor*,  if  P  and  Q  are 
two  local  points  on  a  vortex  line,  then 


(  u  •  7)u  =  I  w|  lim  (11) 

PQ-^oPQ 

where  is  the  variation  u  over  the  distance  PQ. 
That  component  of  ju  which  is  aligned  with  u  acts  to 
stretch  or  compress  the  vorticity  while  the  component 
perpendicular  to  u  causes  the  vortex  line  to  undergo 
instantaneous  rigid  body  rotation.  For  example  if  we 
assume  the  vorticity  is  aligned  in  the  £3  direction  then 


u/3 


du\ 

6x3 


is  the  rotatioh'of  the  vorticity  vector  from  the  ^3  com¬ 
ponent  to  the  wi.  A  similar  interpretation  holds  for 
In  the  remaining  direction  the  vorticity  and 
the  velocity  gradient  are  aligned  and  stretching  or 
compressing  of  the  vorticity  occurs.  This  is  shown 
in  Fig.  2.  The  second  term  on  the  right  hand  side  of 
equation  (10)  is  the  convective  part  of  the  substantial 
derivative. 

Terms  (1),  (2),  and  (3)  of  Eq.  8  can  be  identi¬ 
fied  with  the  convective  part  of  the  nonlinear  terms: 
(1)  and  (2)  can  be  considered  as  transport  terms  and 
(3)  being  a  gradient  production  term.  For  the  ge¬ 
ometry  under  study  term  (1),  transport  by  the  mean 
velocity  is  zero.  Terms  (4),  (5)  and  (6)  are  stretch¬ 
ing/compression/rotation  terms.  Term  (4)  accounts 
for  the  production  of  wiwi  due  to  the  rotation  of  wj 
into  u\  by  This  term,  which  is  non-zero  only 
for  «  =  I,  and  j  =  2  represents  the  action  of  the 
mean  velocity  on  the  fluctuating  vorticity.  Term  (5) 
is  the  production  of  enstrophy  due  to  the  stretch¬ 
ing/compression  and  rotation  of  fluctuating  vorticity 
by  the  fluctuating  velocity  field.  As  will  be  shown 
later,  this  is  an  important  term  near  the  free-surface. 
Term  (6),  the  so-called  mixed  production  term  ac¬ 
counts  for  the  production  of  enstrophy  due  to  the 
stretching  and  rotation  of  the  mean  vorticity  field, 
^3,  by  the  fluctuating  velocity  fields,  u<.  This  term 
represents  the  action  of  the  fluctuating  velocity  field 
on  the  mean  vorticity  field.  The  last  two  terms  are  in¬ 
terpreted  as  viscous  diffusion  and  viscous  dissipation 
of  fluctuating  enstrophy. 

8,  The  Enstrophy  Balances 

The  terms  of  Eq.  (8)  have  been  calculated  and 
are  shown  in  Fig.  3.  They  are  normalized  by  the  ini¬ 
tial  free-surface  velocity  and  the  channel  height.  The 
balance  for  the  total  enstrophy  and  its  components 
near  the  free  surface  are  in  Figs.  4  a-d.  Gains  in  en¬ 
strophy  are  indicated  on  the  right  hand  side  of  these 


Fig,  2,  A  schematic  of  vortex  stretching  and  f'/!  I',  ; 
by  the  fluctuating  velocity. 

figures,  The  individual  terms  give  a  better  indi" n, 
of  the  mechanisms  by  which  the  enstrophy  is  g^n'-r 
•aled  and  destroyed.  .Note  that  the  scale  of  the  ih 
scissa  has  been  changed  between  Fig.  4a  and  Fizs 
4b-d. 

•5.1  Balance  Near  the  Wall 

The  majority  of  the  production  and  dissipation 
of  vorticity  fluctuations  occurs  near  the  wall  wher** 
the  mean  velocity  gradient  and  the  mean  spanwi.s^ 
vorticity  are  large.  In  a  very  thin  region  near  the  wall 
the  diffusion  and  dissipation  of  enstrophy  are  in  h.ii- 
ance.  In  the  region  (-0.8  >  xj  >  -0.95)  the  mix- : 
production  term  (6),  the  production  by  the  mean  .> 
locity  gradient  term  (4),  and  the  turbulent  pr''  '.  1 
tion  terms  (-5)  are  the  dominate  sources  for  flucf  .  i’ 
ing  enstrophy.  in  increasing  importance  away  fr  :• 
the  wall.  These  terms  are  the  rotation/strei'-fnnz 
terms  of  the  balance  equations.  The  gains  ar-  i 
most  completely  balanced  by  the  viscous  dissip-ai.  • 
The  transport/convection  related  terms  are  all  sn, , 
With  the  exception  of  near  the  free  surface,  the  t-rr-  • 
of  the  individual  balances  decrease  wii|, 
distance  from  the  wall.  Near  the  free-siirf,\i-e  th-r- 
a  thin  region  (rj  >  -O.IO)  of  significant  cti.anz-  I 
tails  of  the  individual  balances  near  the  w,all  •  ui 
found  in  Leighton,  et  al.'. 

5.2  Balance  Near  the  Surface 

As  can  be  seen  in  Fig.  4a,  the  viscous  disMi  t 
tion  of  enstrophy  (8)  and  the  turbulent  produ- ^ 

(5)  dominate  in  the  region  -0.10  >  xj  >  -0  3.3  ' 
this  region  there  is  still  some  production  of  flun  ixi 
ing  enstrophy  by  the  mean  velocity  gradient  '>  ■ 
the  mean  velocity  gradient  must  approach  zero  .m  : 
shear  free  surface,  the  enstrophy  production  d  ;•  • 
the  mean  gradient  must  go  to  zero.  Turbulent  .  • 
duciion  of  enstrophy,  the  dominate  production  1; 
anism  for  the  upper  1/3  of  the  channel  decays  rtj 
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for  xj  >'  -0.05.  VUcoui  diffusion  and  turbulent 
transport  act  to  increase  the  enstrophy  in  the  region 
X]  >  -0.10.  Enstrophy  dissipation  balances  the  pro* 
duction  and  diffusion  near  the  top  boundary.  In  the 
thin  region  y  >  -0.08  the  significant  changes  in  the 
terms  of  Cq.  8  can  be  understood  by  examining  the 
component  terms. 

Figs.  4b-d  show  the  level  of  gain  or  loss  for  each 
component  of  enstrophy  near  the  surface.  In  the  up¬ 
per  1/3  of  the  channel  the  transport  term  becomes  rel- 
.Ttively  more  important,  albeit  the  term  is  .still  small 
The  fluctuating  normal  velocity  field  transports  -.i-ii 
and  -,'3u'3  towards  the  surface.  As  the  gradients  m 
and  *3  increase  near  the  surface,  ther'*  is  a  cor¬ 
responding  increase  in  the  enstrophy  transport  (si*.* 
Fig.  1).  The  action  of  the  transport  term  of 
while  always  small,  is  to  reduce  the  level  of  norrtial 
vorticity  near  the  surface. 

Figs.  4b  and  4c  show  that  the  mi.x-d  production 
terms  do  play  a  role  in  the  top  of  the  chanii'd.  but  ap¬ 
proach  zero  as  the  shear-free  surface  is  afiproachei] 
It  is  interesting  to  note  that  the  mi.xi'd  production 
term  for  the  balance  equation  is  negative  in  tlie 
region  (0.0  >  y  >  -0.33)  and  almost  compdetely  bal¬ 
anced  by  the  production  due  to  tlie  mean  velocity 
gradient.  In  the  case  of  wiwi ,  this  term  expres.ses  the 
effect  of  the  rotation  of  the  mean  vorticiiy  vector  by 
the  fluctuating  velocity  field  to  cancel  the  fluctuating 
a.xial  vorticiiy.  This  is  done  at  nearly  the  same  rate 
at  which  vorticity  is  being  rotated  out  of  tlie  mean 
vorticity  into  the  fluctuating  normal  vorticity  by  the 
fluctuating  normal  velocity. 

The  dissipation  and  diffusion  curves  in  Figs  -Ib- 
d  indicates  a  strong  dependence  of  the  diffusion  and 
dissipation  process  on  the  orientation  of  the  vorticity. 
Viscous  diffusion,  which  has  been  negligible  in  the  re¬ 
gion  -O.l  >  y  >  -0.85  becomes  the  dominate  source 
for  the  horizontal  components  of  enstrophy  (wiu/i  and 
.J3W3)  near  the  surface,  (y  >  -0.05).  Due  to  the  re¬ 
flective  nature  of  the  shear  free  boundary,  the  horizon¬ 
tal  components  of  vorticity  go  to  zero  at  the  bound¬ 
ary  and  any  such  vorticity  near  the  boundary  will  be 
influenced  by  the  anti-parallel  image  vorticity  The 
effect  of  the  image  vorticity  or  equivalently  the  shear- 
free  boundary  is  to  increase  diffusion  of  the  vortic- 
ity  in  the  fluid  towards  its  image,  and  m  the  process 
increase  the  dissipation  of  enstrophy.  For  the  nor¬ 
mal  component  (wjwj)  the  process  of  dissipation  and 
diffusion  are  entirely  different.  In  this  case  the  im¬ 
age  vorticity  is  not  anti-parallel  vorticity  promoting 
diffusion,  but  simply  a  continuation  of  the  physical 
vorticity  across  the  shear-free  interface,  The  dissipa¬ 
tion  of  the  mean  square  fluctuating  normal  vorticity 
in  Fig,  4c  is  related  to  the  radial  diffusion  of  normal 
vorticity. 


Fig  3.  The  terms  of  Eq.  8. 


,\lso  apparent  in  Figs.  4b-d  are  the  differences 
in  the  characteristics  of  the  turbulent  production  and 
the  transport  of  enstrophy  between  the  cases  of  the 
horizontal  vorticity  and  the  vertical  vorticity.  With 
the  e.xception  of  the  mean  and  mixed  production 
terms  (4)  and  (6)  the  horizontal  terms  are  similar  In 
the  upper  1/3  of  the  channel,  the  turbulent  prodm' 
lion  of  enstrophy  contributes  positively,  exhibiting  i 
generally  slow  decline  with  increasing  distance  from 
the  solid  wall.  However  in  the  layer  —0.05  >  y  >  -d  I 
the  decline  stops,  and  at  least  for  the  transverse  com 
ponent  of  enstrophy,  there  is  a  slight  increase  in  pro¬ 
duction.  Since  Di  and  fls  approach  zero  at  the  she.ir 
free  boundary,  the  production  of  these  components  •  i 
enstrophy  must  also  approach  zero.  The  n-  u-zero  tur 
bulenl  production  of  the  normal  compoiieni  of  ensir-- 
phy  ,at  the  surface  ,a  result  of  the  shear-free  boun  : 
ary  and  a  normal  strain  field,  which  will  be  discussi- ; 
later. 

The  production  term  (5)  is  further  separated  int 
the  various  stretching  and  rotation  components  1,1 
Fig,  5.  For  all  three  enstrophy  components  the  pr 
duction  due  to  stretching  dominates.  For  the  sirei-  ti 
ing  terms  to  be  relevant,  there  must  be  some  strain 
aligned  with  the  vorticity.  While  this  strain  is  not  .»■ 
counted  for  by  the  existence  of  the  anti-parallel  imag- 
vortex,  numerical  experiments  involving  interacting 
vortex  rings  or  vortex  tubes  indicate  a  rapid  incr**  t- 
in  the  component  of  strain  aligned  with  the  vorti  •• 
during  the  vortex  reconnection.  A  similar 


Fig.  4a.  Enstrophy  balance  in  a  lurbul-'tu  '.pf'ii 
channel  flow  near  the  free  surface. 


Fig.  4c.  .Mean  square  fluctuating  normal  vorticity 
balance  in  a  turbulent  open  channel  flow  near  the  free 
surface. 


Fig.  4b.  Mean  square  fluctuating  axial  vorticity  bal¬ 
ance  in  a  turbulent  open  channel  flow  near  the  free 
surface. 

process  may  be  relevent  when  the  horizontal  compo¬ 
nents  of  vorticity  approach  the  top  surface  The  flow 
field  transporting  the  vorticity  to  the  free  surface  may 
also  be  responsible  for  the  strain. 

Although  the  rotation  terms  are  small  and  de¬ 
clining  as  the  surface  is  approached,  the  terms  repre¬ 
senting  the  rotation  of  vorticity  from  the  axial  to  the 
transverse  direction  and  from  the  transverse  to  the 
.axial  have  small  increases  near  the  boundary.  This 
may  be  due  to  the  rotation  of  small  eddies,  having 
vorticity  components  parallel  to  the  surface  by  large 
eddies. 


Fig.  4d.  Mean  square  fluctuating  transverse*  vortiriis 
balance  in  a  turbulent  open  channel-  flow  n**ar  the  fr**** 
surface. 

The  enstrophy  balance  in  the  vertical  (rj)  direction 
shows  the  effect  of  attachment  of  the  normal  vortn- 
ity  to  the  free  surface.  The  turbulent  production  -  i 
wjwj  varies  slightly  in  the  upper  1/3  of  the  channel 
However,  as  can  be  seen  in  Fig.  5  the  production 
mechanism  changes  with  increasing  distance  from  the 
shear-free  boundary.  Near  the  surface,  xj  >  -0  0" 
the  stretching  of  the  vorticity  is  important.  Below 
(xj  <  -0.08)  the  gain  in  enstrophy  is  due  to  rotation 
of  the  fluctuating  vorticity  field,  and  stretching  hv« 
become  negligible.  The  source  of  the  strain  near  tlie 
free-surface  is  probably  the  flow  induced  by  the  !•  w 
pressure  within  the  vorticity  attached  to  the  surf.a-- 
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From  figure  5,  it  it  apparent  that  the  term  in 
Eq.  8  meet  reepotuible  for  the  production  of  en- 
strophy  i«  traniverae  etretching  of  traniverie  vortie- 
ity,  W3W3  Hi  Fipit.  6  U  a  piot  of  thia  property  for 
the  piane  Z2  s  >0.06.  The  turbuient  itructurea  or 
eventa  reaponaible  for  generating  thia  component  of 
the  enatrophy  are  localized  and  intenae.  The  contour 

intervale  are  10w3W3|^.  To  more  clearly  identify  the 
structures  associated  with  the  localized  production, 
histograms  of  the  level  of  production  by  turbulent 
stretching  have  been  determined  and  are  contained 
in  Leighton^,  et.al.  From  the  histograms  it  was  deter* 
mine  that,  at  this  plane  (zj  =  >0.05)  that  approxi¬ 
mately  3  -  8%  of  the  area  was  responsible  for  30  -  60% 
of  the  turbulent  enatrophy  production.  This  behav¬ 
ior  is  typical  for  for  all  nine  components  of  term  (.5). 
Eq.  8.  Thia  conclusion  was  verified  by  calculating 
the  skewness  and  flatness  for  these  production  terms. 
For  each  of  the  nine  cases  the  skewness  and  flatness 
factors  were  large  (10-  10*  for  skew  and  10*  -  10^  for 
flatness). 

Based  on  the  observation  of  production  events  lo¬ 
cal  in  space,  and  the  behavior  of  the  average  fluctuat¬ 
ing  vorticity,  two  models  of  the  interaction  of  vorticity 
with  the  surface  were  developed. 


6.1  Two  models  for  the  production  events 

The  two  proposed  models  of  the  interaction  of 
vorticity  with  the  free  surface  can  ,  be  described  as 
‘the  spin  model’  and  ‘the  splat  model'.  These  mod¬ 
els  are  sirnilar  to  and  have  been  infleuenced  by  the 
models  presented  by  Koh  and  Bradshaw^  and  Hunt*. 
Although  the  names  are  meant  to  be  suggestive  of 
the  kinematics  of  the  enatrophy  producing  events,  the 
models  are  similar  to  the  models  by  the  same  name 
for  the  pressure  source  terms  described  in  Ref,  9. 

In  the  'spin'  model,  a  vortical  structure  originat¬ 
ing  in  the  buffer  layer,  tilts  downstream  and  is  at¬ 
tached  to  the  free  surface  in  a  quasi-stable  configu¬ 
ration.  This  model  is  relevant  when  discussing  the 
vertical  component  of  enatrophy.  These  structures, 
which  can  be  liken  to  vortex  tubes  connected  to  the 
free  surface  have  been  observed  in  flow  visualizations 
of  the  DNS  dataset  of  the  open  channel  flow.  Within 
these  vortex  tubes  there  is  a  significant  'down  draft’. 
It  is  conjectured  that  this  down  draft  is  responsible 
for  the  transport  of  enatrophy  away  from  the  surface 
and  provides  the  strain  necessary  to  maintain  these 
structures.  Were  it  not  for  these  down  drafts,  and 
the  vorticity  intensification  they  provide,  these  vorti¬ 
cal  structures  would  eventually  diffuse  and  cancel. 


Fig.  6.  Turbulent  production  of  fluctuationg  en»ir'^>- 
phy  by  stretching  (  heavy  lines)  and  rotation  (light 
lines),— 

The  eventual  diffusion  has  been  seen  in  a  simulation 
of  vortex  reconnection  in  which  there  was  no  down 
draft.'®  In  this  case  there  was  radial  diffusion  of  vor¬ 
ticity  from  the  core  near  the  intersection  of  the  vortex 
ring  and  the  shear-free  surface. 

In  the  ‘splat’  model  a  volume  of  fluid  ‘contai!. 
ing’  vorticity  impacts  the  free  surface,  Ie.ading  ’ 
stretching,  intensification  and  cancellation  of  vori. 
ity.  The  fluctuating  normal  velocity  field  proni' 
the  increase  of  horizontal  component  enstrophy  n--ir 
the  free  surface  by  two  mechanisms;  transport  .in  ; 
vortex  stretching.  Due  to  the  large  gradient  in  - 
and  near  the  top  boundary,  the  vertical  velo-  :■ . 
fluctuations  provide  an  efficient  means  of  transporti;  i 
vorticity  to  the  free  surface.  As  the  turbulent  •‘■i  l. 
approaches  the  top  boundary  the  vertic.ii  (n'TU 
is  deflected  into  the  horizontal  direction  During  'H  - 
’splat'  event,  the  vorticity  contained  in  a  lurbiil-  !  ■ 
eddy  interacting  with  the  top  surface  will  be  roti!-  • 
parallel  to  the  boundary  and  stretched.  The  vr.r*. 
ity  in  the  fluid  near  the  free  surface  will  be  subj**'  ••  • 
to  a  large  strain  and  will  have  a  large  gradient  ]>.■ 
to  the  proximity  of  its  reflected  image.  This  leads  • 
increased  production  of  fluctuating  enstrophy  by  ti  • 
fluctuating  velocity  fields  and  to  large  enstrophy  I.' 
sipation  or  vorticity  cancellation.  Other  details  ah-  .- 
this  type  of  model  can  be  found  in  Ref.  9. 
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Cfladitioniliampllnc  and  EniamHU  Avarayin| 

To  verify  tbtt  the  two  modeli  tre  quelitetiveiy 
correct,  the  open  channel  flow  waa  conditionally  iam* 
pled  on  the  plane  n  »  -0.06.  The  condition  need 
wai  that  the  local  value  of  the  enetrophy  production 
by  vortex  itretching  ( rotation  wai  not  coneidered  ) 
exceeded  the  mean  value  by  a  factor  of  ten,  and  that 
the  production  wee  a  local  maximum.  The  enetro¬ 
phy  production  at  the  locations  determined  by  the 
first  condition  accounted  for  approximatedly  50%  of 
the  production  by  stretching  and  rotation,  but  only 
about  2%  of  the  area.  For  the  spanwise  component. 

only  detected  events  with  negative  vorticity 
were  considered.  There  were  approximately  an  or¬ 
der  of  magnitude  more  events  with  negative  vorticity 
contributing  to  the  turbulent  production  of  spanwise 
enstrophy  than  events  with  positive  vorticity.  For  the 
remaining  components,  only  events  with  positive  vor¬ 
ticity  were  considered.  In  this  case  there  were  as  many 
positive  as  negative  vortical  events  detected. 

Certain  properties  associated  with  these  events 
have  been  ensemble  averaged  and  are  shown  in  figures 
7  -  10.  The  properties  are  averaged  on  a  vertical  line 
through  the  detected  events.  The  ensemble  averaged 
vorticity,  shown  in  Fig,  7  are  larger  than  those  shown 
in  Fig.  1,  almost  by  a  factor  of  ten,  an  indication  that 
these  events  are  substantial  contributors  to  the  rms 
vorticity. 

7.1  The.  Spin  Event 

As  discussed  above  the  spin  model  involves  a  vor¬ 
tical  structure,  like  a  vortex  tube  attached  to  the  free 


surftes  Md  MUndiag  into  tb«  buiCsr  Uytr,  The  «x- 
istMce  of  subftMtial  «OMmblo*avtr»f«d  vertical  vor¬ 
ticity  in  regione  far  from  the  free  surface,  a«  seen  in 
Fig,  7  supports  the  conceptual  model  of  an  attached 
vortex  tube. 

Fig,  8  is  the  ensemble-averaged  vertical  velocity. 
In  the  ease  of  the  normal  component  of  enstrophy, 
the  negative  velocity  will  stretch  the  attached  vertical 
vorticity.  This  is  the  down-draft  within  the  vortex 
tube  mentioned  earlier.  The  peak  negative  veiocity 
occurs  at  *j  -0.15.  The  strain  rates  ^  are  shown 
in  Fig,  9.  The  peak  in  the  normal  strain  occurs  at 
the  surface,  again  consistant  with  the  proximity  of  (he 
vortex  stretching  with  the  free  surface.  The  strain 
and  the  resulting  vortex  stretching  are  negative  for 
rj  <  -0.15.  This  may  be  due  in  part  to  the  tilting  ot 
the  vortical  structure  in  the  main  part  of  the  flow 

The  final  figure  is  the  rate  of  dissipation  of  ''n- 
atrophy  for  each  of  the  components.  The  normal  v.  r 
ticity  is  largely  unaffected  by  its  image  and  has  a  low 
level  of  dissipation.  Due  to  this  small  dissipation,  it 
is  conjectured  that  these  events  will  have  a  longer  du¬ 
ration  than  the  ‘splat'  events. 

7.2  Ihe_S2|ft.t.,Ev?Jl^ 

Although  the  conditionally  sampled  results  ar** 
consistant  with  both  the  model  and  the  RMS  vorticitv 
distributions,  they  donot  explain  the  greater  intensity 
of  the  W3  components.  The  location  of  the  maximum 
transverse  vorticity  is  above  the  conditioning  plan 
y  sr  -0,05  while  the  maximum  for  the  sxial  vorticity 
is  below,  This  characteristic  and  the  amplitude  of  ih- 
maximum  are  consistant  with  their  rms  values 


Fig,  6,  Enstrophy  production  due  to  transverse  stretching  of  transverse  vorticity  (  wswslj**).  Planview  of 
plane  y  =s  -0,06.  Contours  are  for  s  lOujus^. 
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Fig,  8  ii  t¥«  vertical  velocity  auociated  with  thoe 
event*.  For  the  rieinf  turbulent  eddy  of  the  model 
the  vertical  velocity  is  poeitive,  reeulting  In  a  'iplat* 
at  the  surface.  The  strain  rates  are  shown  in 
figure  9.  The  transverse  strain  is  larger  is  the  axial 
strain.  This  is  consistent  with  the  larger  transverse 
turbulent  production  and  vorticity. 

Fig.  10  is  the  rate  of  dissipation  of  enstrophy  for 
each  of  the  components.  The  terms  in  xi  and  xs  direc¬ 
tion  are  similar.  The  dissipation  occurs  very  close  to 
the  boundary  and  is  very  intense:  approximately  an 
order  of  magnitude  larger  than  the  mean  value.  Due 
to  the  reflective  nature  shear-free  boundary  there  ex¬ 
ists  an  image  vorticity  across  the  boundary  resulting 
in  annihilation  of  vorticity  near  the  shear-free  bound¬ 
ary. 

8.  Conclusions 


The  enstrophy  balance  equations  have  been  eval¬ 
uated  for  a  direct  numerical  simulation  of  turbulent 
open  channel  flow,  The  free  surface  has  been  mod¬ 
eled  as  a  shear-free  boundary.  Near  the  free  surface 
there  are  two  dominate  sources  of  enstrophy;  Trans¬ 
port  due  to  velocity  fluctuations  and  production  due 
to  the  vortex  stretching  and  rotation  by  the  velocity 
fluctuations.  At  the  free  surface,  diffusion  is  the  dom¬ 
inate  source  of  enstrophy.  The  gains  in  enstrophy  are 
balanced  by  enstrophy  dissipation.  When  the  compo¬ 
nents  of  the  total  enstrophy  are  considered  individu- 
ally,  the  components  parallel  to  the  surface,  (-iwi  and 
wiwi)  are  found  to  be  responsible  for  the  variation  in 
the  balances  near  the  surface.  The  balance  equations 
for  wjwj  show  little  variation  near  the  surface. 

Two  models,  referred  to  as  the  'spin'  and  'splat' 


models  can  be  used  to  explain  this  behavior.  The 
splat  model  Is  similar  to  the  model  discussed  in  Hunt* 
in  which  an  eddy  Is  swept  to  the  surface.  The  result¬ 
ing  redistribution  of  normal  momentum  into  horizon¬ 
tal  momentum  will  stretch  and  intensify  any  vorticity 
transported  in  the  eddy.  Due  to  the  image  vorticity 
resulting  from  the  reflective  nature  of  the  free  surface, 
the  intense  vorticity  will  quickly  cancel  with  its  image 
and  the  enstrophy  dissipation  will  be  enhanced. 
Fundamental  to  the  ‘spin’  model  is  a  stable  vor¬ 
tical  structure  originating  in  the  buffer  layer  and  ex¬ 
tending  to  the  free  surface.  The  mild  variation  of 
the  normal  balance  is  a  result  of  the  quasi-stable  n<a- 
lure  of  the  attached  structure.  Without  close  ,inii- 
par-allel  image  vorticity,  these  structures  do  not  dis¬ 
sipate  quickly.  The  downflow  within  these  vortic.al 
structures  result  in  sufficient  strain  to  maintain  th*- 
vorticity  against  the  effects  of  viscous  diffusion 

A  large  portion  of  the  production  of  enstrophv 
near  the  free  surface  is  due  to  these  highly  localiz.-d 
intense  events,  The  area  responsible  for  this  produc¬ 
tion  is  very  small.  Approximately  3  -  6%  of  the  area 
on  the  plane  xj  =  -0.05  is  responsible  for  30  -  WA 
of  the  production.  Conditionally  sampling  and  en¬ 
semble  averaging  the  production  events  verified  that 
their  characteristics  are  consistant  with  the  propose.) 
models. 
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Fig,  7.  Ensemble  averated' vorticity,  conditioned  by  Fig,  8,  Ensemble  averaged  verUcal  velocity, 

the  requirement  that  >  lOwtWi 
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Fig.  9.  Ensemble  averaged  strain. 
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